Discrete Green's functions for products of regular graphs 

Robert B. Ellis* 

Department of Mathematics, Texas A&M University, 
College Station, TX 77843-3368, USA 



Abstract 

Discrete Green's functions are the inverses or pseudo-inverses of combinatorial Laplacians. 
We present compact formulas for discrete Green's functions, in terms of the eigensystems of 
corresponding Laplacians, for products of regular graphs with or without boundary. Explicit 
formulas are derived for the cycle, torus, and 3-dimensional torus, as is an inductive formula 
for the t-dimensional torus with n vertices, from which the Green's function can be completely 
determined in time 0{t-n?~^^^ logn). These Green's functions may be used in conjunction with 
diffusion-like problems on graphs such as electric potential, random walks, and chip-firing games 
or other balancing games. 
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1 Introduction 



Discrete Green's functions arc the inverses or pseudo-inverses of combinatorial Laplacians, which 
govern diffusion-like problems on graphs such as random walks, electric potential and chip-firing 
games or other balancing games. Just as Green's functions in the continuous case depend on the 
domain and boundary conditions, discrete Green's functions are associated with the underlying 
graph and boundary conditions, if any. Just as in the continuous case, a new set of discrete Green's 
functions must be determined for each new class of graphs. Certainly, a discrete Green's function 
can be determined by brute force (pseudo-) inversion of the corresponding Laplacian, but this is no 
advancement toward compact or closed- form functions. 

In this paper, we develop such compact formulas for discrete Green's functions on products of 
simple graphs with or without boundary. Section 2 presents the necessary definitions and back- 
ground. Section 3 illustrates these definitions by deriving the Green's function for the cycle. In 
Section 4, we derive formulas for products of regular graphs with or without boundary. In Section 5, 
we illustrate the case of products of regular graphs without boundary by deriving the Green's func- 
tion for the t-dimensional torus and addressing its computational complexity, with explicit Green's 
functions given for t = 2 and < = 3. Finally, the relationship between Green's functions and the 
hitting time of a random walk is explained and illustrated in the case of the 2-torus in Section 6 

* Research supported in part by NSF grant DMS-9977354. Email address: rellis@math.tamu.edu. 
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2 Preliminaries 



The basic definitions follow those of [4]. Let F = (V, E) be a simple connected graph. Let x,y € V 
be arbitrary vertices. The adjacency matrix A is defined by A(x, y) = xi^ ~ y). The diagonal degree 
matrix D is defined by D{x, y) = x{x ~ y) ■ dx, where dx is the degree of x in F. The volume of a 
graph is vol{T) ~ X^uei/ The Laplacian, i = Z? — A, of F is 



L{x,y) 



dx, 
-1, 
0, 



if X = y 
if X ^ y 
otherwise. 



The normalized Laplacian, C — D ^/^LD is 



1, 



C{x,y) = 
The discrete Laplace operator A is 

A(x,y) 



-l/y/dxdy, 
0, 



ii X = y 
ii X y 
otherwise. 




if a; = y 
if a; ~ y 
otherwise. 



In general, the following relations hold between L, £, and A: 

C = D-^'^LD-^'^ 

L = DA. 



For S* C we define the Dirichlet versions of Ls, Cs, and A5 as the results of deleting the 
rows and columns corresponding to y \ S* from L, C, and A, respectively. Without loss of generality, 
we assume that both F and the subgraph generated by S to be connected. By abuse of notation, 
we use S to refer to the subgraph generated by in F. Additionally, we say that the orthonormal 
eigensystem of S is the orthonormal eigensystem {(Aj, (f)j) : j S J} of the real symmetric matrix Cs', 
where J = {0, . . . , |F| - 1} with Aq = when 5* = V', and J = {1, . . . , \S\} with Ai > when S j^V. 
In either case, we order the eigenvalues by Xj < Aj+i for all subscripts in range; basic properties 
such as < Aj < 2 for all j are summarized in [3]. In particular, the orthonormal eigenvectors can 
be chosen to have real entries. 

When S CV, As, Ls and £5 arc invcrtiblc, and the Green's function G and normalized Green's 
function Q are determined by the relations 

AsG = GAs = Ls and 

CsQ = gCs = Ls. (1) 

We can tie these relations to random walks as follows. Let P = [P{x, y)] be the transition probability 
matrix for the simple irreducible transient random walk on S with absorbing states V\S, where the 
probability pxy of moving to state y from state x is \/dx if x and y are adjacent and otherwise. 
Then As = L - P, and (/ - P)-^ = I + P + P'^ + ■ ■ ■ gives 

G(x,y) = ^P„(rE,y), (2) 

n 



2 



where Pn{x,y) is the n-step transition probabiHty matrix (cf. [11, p. 31]). Sec [1. 7] for definitions 
and results on random walks. 

On the other hand, when S = V, the sum in (2) does not converge, and we require an alternate 
definition of the Green's function. Since C is not invertible, its corresponding normalized Green's 
function Q is defined by 

\j>0 ^ 

This definition of Q is equivalent to the two relations 

gC = eg = I - Pq = I ~ 4>o(l3*o and 

GPo = 0, (4) 

where Pq = 0o</'o projection of the orthonormal eigenvector 4>q corresponding to eigenvalue 

0. Since LI = 0, we have 0o = -D^/^l/HZJ^/^lH, and so (jjoix) = ^ydx/vol{T). It is also important 
to note that when g is invertible, the definitions in (1) and (3) are equivalent. Furthermore, g is 
related to the so-called fundamental matrix 

oo 

Z{x,y) ■.= ^{Pn{x,y) - TTy), (5) 

n=0 

where tt is the stationary distribution of the random walk on F, by the equation g = D^^^ZD^^^^; 
simply verify that D^l^ ZD~^I^ satisfies (3) and (4). See Chapter 3, p. 17 of [1] for relationships 
between Z and hitting times. A more complete exposition on discrete Green's functions in the 
context of spectral graph theory appears in [8] . 



3 Green's function for the cycle 

In this section wc illustrate the preceding definitions in the case of the cycle, which has no boundary. 
The techniques developed here will be used in later sections to construct the Green's functions for 
higher dimensional tori. 

By the cycle Cm, we mean the 2-regular connected graph with vertices {0, 1, . . . , m — 1}, where 
m > 3. The various Laplacians are related by A = £ = L/2. Applying the definition in (4), the 
normalized Green's function g satisfies 

gC = eg = I — —J and 

m 

gj = 0, (6) 

where (f>o{x) = i/l/m, and J is the m x m matrix of I's. Because the cycle is vertex-transitive, the 
values C{x, y) and g{x, y) depend only on the distance \y — x\ between x and y, and the following 
definition of g{a) is well-defined: 

g{a) := g{x,y), ifa=|y-a;|. (7) 
We are ready to derive the Green's function for the cycle. 
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Theorem 1. Let m > 3. For < x,y < m ^ 1, the cycle Cm has normalized Green's function 

g{x,y)=^- £^ l-\y^x\+^-^ (8) 

bm m 

Proof. From (6) and (7), wc have the reeurrenee 

2g{x, y) - g{x, y-1)- g{x, y + 1) = ( ^ 1 = ^ or 



2g{a) -g{a-l) -gia+1) = 



-2/m, X =/= y, 

2-2/m, a^O, 
-2/m, a > 0, 



provided that we define g{—l) = 5(1) for shiiphcity of representing the case a = 0. The foUowing 
recurrence on differences results: 

g{a + l)-g{a) = g{a)-g{a^l) + - - 2xia = 0). (9) 

m 

The second constraint in (6) determines that the sum of g across any row must be 0; i.e., 

771— 1 

a=0 

We now solve the recurrence on differences, starting with g[l) — g(0) hy setting a = in (9). 

i (5(1)- 2(7(0) + = 5(1) -5(0) 

= 1-1. (11) 
m 

Resolving (9) with base case given by (11), we have 

5(a + l)-5(a) = 5(a)-5(a-l) + - 

m 

= 5(a-l)-5(a-2) + 2- — 

m 

m 

^ 2 

= l + a-—. (12) 

m m 

Having derived a simple recurrence from the recurrence on differences, we proceed to determine 
5(a). Resolving (12) yields 

5(a) = 5(a-l) + l-l + (a-l) — 
m m 

= 5(a-2) + — -2+[(a-l) + (a-2)] — 
TO m 
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a L 



e(0) + --a+-y 

m m ^ — ^ 



m TO 

k=0 



= e(0)-a+— . (13) 

TO 

Now applying the row sum constraint (10) allows us to compute the value of ^(0). 



m-g{0) = J2 



m — 1 / 

[a- 

a=0 ^ 



a^ 

TO 



TO 



^(0) = (!!l±il(!Ii^ . (14) 

DTO 

Plugging (14) into (13) and letting a—\y~x\ achieves the desired result. □ 



Because (3) also gives Q for Cm, we have a whole class of identities formed by choosing any 
orthonormal eigenbasis for Cm and equating (3) with (8). We give one such well-known basis now, 
which arises naturally from the consideration of circulant matrices (cf. [6]). 

Lemma 1. For to > 3, define (\)j{x) := cxp ■ Then {(l — cos [271] /m), 4>j) ■ < j < to} 

is an orthonormal eigensystem of C for Cm ■ 



The proof is omitted, but follows from a straightforward verification of the necessary conditions. 
Theorem 2 follows by combining Theorem 1 with (3) using the orthonormal eigenbasis of Lemma 1. 

Theorem 2. Let to > 3 and let < x,y < m. Then 

J_ exp {{2Trij/m){y - x)) _ (to + 1)(to - 1) _ . _ _^ iv " 
TO ^-^ 1 — cos {2Trj /m) 6to to 



4 Green's functions for products of regular graphs 



Theorems 4-5 of [3] give a contour integral formula for the Green's function of the Cartesian product 
of two regular graphs with boundary, provided a certain generalized Green's function is known for 
each factor graph. In this section, we extend these results to include the cases where one or both 
of the factor graphs is without boundary, and provide simplified working formulas requiring the 
generalized Green's function of one graph and the eigensystem of the other. The original technique 
of [3] recovers the Green's function of the product graph as the residues of a certain contour integral 
whose contour contains as poles all of the eigenvalues of one of the graphs, and none of the negatives 
of the eigenvalues of the other graph. When both graphs are without boundary, both have an 
eigenvalue of 0, and the resulting order 2 pole requires a modified contour integral formula. 

Recalling the definitions in Section 2, Let F = {V,E) and F' = {V',E') be regular graphs of 
degree d and d' with specified vertex subsets S and S" and normalized Dirichlet Laplacians Cs and 
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C'gi, respectively. For any a G C, define Qa to be the symmetric matrix satisfying tlie relation 
(£5 + a)Qa = /s, if S' C and the relations 



'a 



Is - Po, 



and 



GaPo = 0, 



(15) 



if S ~ V {Cs is singular). In cither case, this is equivalent to 



Gaix,y)= p— 0j(x)0j(y), 

Xj>0 ^ 



(16) 



where the 0j's are the orthonormal eigenfunctions of Cs associated with the eigenvalues Xj. In 
particular, Qa is a rational function of a. The analogous definitions of Q'^, <t>ij and are made for 

r. 

The choice of S and S' induces the specified vertex set S' x 5" in the Cartesian product F x F'. 
This product of F = {V, E) and F' = [V , E') has vertex set {{v, v') : v e V, v' € V'} and edges of the 
form {{v,v'), {v,u')} or {{v,v'), {u,v')} where {u,v} G E, {u',v'} G E' . Let be the normahzed 
Laplacian of F x F', with corresponding normalized Dirichlet Laplacian J~-g-^gi ■ We abuse notation by 
referring to S, S' and 5 x 5" as the graphs they induce. If £3 and Cg, have orthonormal eigensystems 
{{Xj,4ij) : j G J} and {(AJ., 05c) : k G K}, respectively, then when S and S' are regular of the same 
degree, ^sxS' orthonormal eigensystem 



where $j.fc(f,f') := 4'j{v)(j}'f,{v'). We begin with the case in which S has boundary. 
4.1 At least one graph has boundary 

Suppose S C.V , so that S generates a connected subgraph with boundary in F. Wc allow S' C V' 
to generate an arbitrary connected subgraph in F'. We consider two cases, where the factor graphs 
are regular of the same degree or regular of different degrees. 

First, suppose S and S' be regular of the same degree, and let C denote a simple closed contour 
in the complex plane consisting of all a G C satisfying |q; — 1| = l + Ai/2. The contour C is designed 
to enclose all of the A'j,'s and none of the — Aj's. We have the following generalization of Theorem 4 



Theorem 3. Let S and S" be induced subgraphs of T = (V, E) and F' = {V' , E'), respectively, which 
are both regular of degree d. Let S V and S' C V' . The normalized Green's function G of the 
Cartesian product S x S' with Dirichlet boundary condition is 



Proof. Combining the product graph eigensystem in (17) with the formal definition of G in (3), we 
have 




(17) 



of [3]. 
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\s\ 



3 = 1 k 



(A,+a)(A',-a) 



ga{x,y)g'_a{x',y)da, by (16). 



□ 



Note that the above contour integral picks up exactly the residues at a = A'j,. For example, the 
residue at a = AJ, is exactly 

^ ^ J^MfMl^lpZI ^ J2 , ^'K(^')M)-o.My)- (18) 

For convenience, we assign the term of the residue in (18) corresponding to K to A'^^ . This observation 
gives us the computational formula for G in the following corollary, which in practice may be applied 
to yield a closed formula for a specific product graph, or to generate in conjunction with (3) a non- 
trivial identity involving Qa and arbitrary orthonormal eigensystems of Cs, Cs' and -Cg^^, . 

Corollary 1. Under the same conditions as in Theorem 3, we have 

G{{x, x'l (y, y')) = 2 ^ </>L.(a:') WI^a;. {x, y). 

k 



Second, suppose we have the same conditions as Theorem 3, except S and S' are regular of degrees 
d and d', respectively. The normalized Dirichlct Laplacian C^y,g, has orthonormal cigcnsystcm 

{(^A, + -^A'„^:.eJ,fceA| (19) 

The poles of Qa/d sre at a = —dXj, and the poles of Q-a/d' ^tre at a = d'Aj.. Let C denote a simple 
closed contour in the complex plane consisting of all a G C satisfying \a — d'\ = rf' + rfAi/2; thus 
C contains all of the <i'A'j.'s but none of the — rfA^ 's. Wc obtain the following minor extension to 
Theorem 5 of [3]. 

Theorem 4. Let S and S' be induced subgraphs ofT^ (^i-^) '^^'^ T' = {V-iE!'), respectively, where 
r is regular of degree d and F' is regular of degree d' . Let S C. V and S' C V' . The normalized 
Green's function G of the Cartesian product S x S' with Dirichlet boundary condition is 

G((x, x), {y, y')) ^ ^^^^ Ga/dix, v)Q'-a/d' i^' ^ y')da. 



Proof. Combining (19) with (3), we have 

G((x,x'),(2/,y')) = dxXd'K '^'^'^''''''^^'''^^'^'^ 

d + d' /■ ^ >p 0, {x)JM<P'k {x'WM 

2^^* Xc^tt {dXj+a){d'\'^-a) 
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d + d' 

2TTidd' 




{X,+a/d)iX,-a/d') 



d + d' 
2TTiddf 



Jc 




□ 



Analogous to Corollary 1, by inspecting the residues of the above contour integral at all values 
a — d'A'j,. we have the following. 

Corollary 2. Under the same conditions as Theorem 4, we have 



4.2 Neither graph has boundary 

Here we consider the case oi S ~ V and S' ^ V , so that S* x 5' is the entire product graph F x F'; 
recall that the normalized Green's function for F x F' is non-invertible. Let m ~ \V\ and n = \V'\, 
and consider first the case in which F and F' have the same degree. 

The orthonormal eigensystems of C and C are {{4>j, Aj) : < j < m — 1} and {((/>fc, AJ,) : < 
fc < n — 1}, respectively. Let C denote a simple closed contour in the complex plane consisting of 
all a G C satisfying |a — (2 + A'^/2)| = 2. The contour C is designed to enclose A'^, . . . , A^_^ and 
none of — Aq, . . . , —^m-i or Aq. 

Theorem 5. Let F and F' be connected regular graphs of degrees d without boundary. With the 
notation above, the normalized Green's function G of the Cartesian product F x F' is 



Proof. The eigenvector (j)o corresponding to eigenvalue is determined by (j)o{x) = ^ {dx/vol{r)). 
Combining the formal definition of G in (3) with (17), and noting that dy = dv> = d, vol{A) ~ d-m, 
and vol{A') — d - n, we have 




k 







(l)j{x)<i)j{y)(j)k{x')(j)k{y') 
(A,+a)(Al,-«) 



da 
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= - [ g^ix,y)g'_^{x',y')da + -g{x,y) + -g\x',y') by (16). 

□ 

Analogous to Corollary 1, inspecting the residues of the above contour integral at X'l, . . . , X'„_i 
yields the following corollary. 

Corollary 3. Under the same conditions as in Theorem 5, we have 

G{{x, x'), (y, y')) = 2 ^ ^U^^O^U^O^a^ v) + Six, y) + -^'(x', y'). 

k=l 

Now suppose wc have the same conditions as Theorem 3, except that the graphs F and V are 
regular of degrees d and d' , respectively, and -Cp^^p, has orthonormal eigensystem given by (19). Let 
C denote a contour in the complex plane consisting of all a S C satisfying \a — {d' + d'X[/2)\ = d' . 
Thus C contains d'X[, ■ ■ ■ , d'X[^_i, but neither — dAg, . . . , ~dXm-i nor cZ'Aq. We obtain the following 
theorem. 

Theorem 6. Let T and T' be connected regular graphs without boundary of degrees d and d' , respec- 
tively. With the notation above, the normalized Green's Junction G of the Cartesian product F x F' 
is 

G{{x, x'l {y, 2/')) = ^ I G./Ax, y)g'.^,a' (^', y')da + ^^(x, y) + '^±^g\x' ,y'). 



Proof. Combining the definition of G in (3) with (19), we have 

d + d' 
dXn + d'X 



G{{x,x'),iy,y')) = J2J2 dx +d'X' '^jM^^^')'^,Ay,y') 

j = l k=l 



ni-1 ^ _|_ ^/ d + d' 

^73~*J'0(^' ^')^xo(y, y') + ^70^*"^^^^' ^')*o,fc(y, y') 



j=i ^ k=i 



d + d' f dd'(t)j{x)4>j{y)4>k{x')4>k{.y') 



2'Kidd' Jc {dXj + a){d'Xk - a) 



I d + d'./d^ d + d'^y 

d + d' f ^ , , ^, , . , d + d! , d + d! ^1, I ,. 



□ 



Analogous to Corollary 1, by inspecting the residues of the contour integral at d'X'i, . . . , d'X'^_i, 
we have the following. 

Corollary 4. Under the same conditions as in Theorem 6, we have 

G((x, x'l (y, y')) = ^ E <l>'kix')-^)g,,yj,{x, y) + ^^g{x, y) + ^^'(x', y'). 
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5 Example Green's functions for products of cycles 



Application of the results in the previous two sections to specific examples is limited only to cases 
in which the necessary raw materials can be computed. The results in Section 4.1 for the product 
where at least one graph has boundary require the eigensystem of one graph and the generalized 
Green's function Qa of the other graph. Although as written, the results assume that the generalized 
Green's function is known for the graph without boundary, and the eigensystem is known for the 
graph without boundary. Theorems 3 and 4 can be easily re-derived for when the reverse is true. 

The results in Section 4.2 for the product of two graphs without boundary require knowledge 
of the eigensystem of one graph, the generalized Green's function of the other graph, and the 
Green's functions Q and Q' of both graphs. Thus computation of the normalized Green's function for 
any graph can be done via any decomposition of the graph into factor graphs where this information 
is known. This observation can be particularly useful in computing the normalized Green's function 
inductively where each additional factor graph is from a specific family. 

5.1 The torus Cm x C„ 

Following Corollary 3, determination of the Green's function of the torus requires the Greens function 
Q and generalized Green's function Qa of the cycle. In obtaining a compact formula for the torus, 
it is critical to simplify Qa as much as possible before incorporating it into GoroUary 3. 

Theorem 7. Let to, n > 3. For < x,y < m — 1 and < a;', j/' < n — 1, the torus Cm x Cn has 
normalized Green's function 



where T and U are the Chebyshev polynomials of the first and second kind, respectively. 

Note that the formula depends only on the distances between y and x and between y' and x' , 
which is expected due to the translational symmetries of the torus. The first step of the proof is 
to obtain a closed form for Qa for the cycle Cm {Q was determined in Theorem 1). The proof of 
Theorem 7 is deferred until after Cor. 5. 

Theorem 8. Let m > 3. For a cycle Cm with vertices 0,1,..., to — 1, complex a ^ 0, and 
< x,y < m ~ 1, the generalized Green's function Qa satisfies 




Tm/2-\y-x\{2 - COS (27rfc/n)) 



(1 — COS (27rfc/n))(3 — cos (27rfc/rt))C/j„/2-i(2 — cos {2'Kk/n)) 




Ga{x,y) 



2 



+ 




TO (r + r 



1 



where 2(1 + a) = r + r 
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Proof. Because Cm is vertex-transitive, Qa{x, y) depends only on the distance mindy — xl, to— — x|) 
between x and y. Therefore define GaiO') ■— Qaix, y), where a — \y ~ x\; this induces the additional 
relation Qa{ci) — Qa{m — a) for 1 < a < to — 1. From (15), we have 

X{x = y)-— = {£s + a)ga{x,y) 

TO 

- i (2(1 + a)g^{x,y) ^g^{x + l,y) - G^ix - 1, y)) 
= i((r + r-i)g„(.T,y)-e„(a- + l,y)-g„(x--l,y)), 

where Cs is the normalized Laplacian of Cm- We can rewrite this as 

2 1 

Ga{x - l,y) - rga{x,y) = 2x{x = y) + -{ga{x,y) - rg^ix + l,y)), 

m r 

which for a > becomes 

2 1 

ga{a + 1) - r g^ia) = h - (^a(a) - r ^/^(a - 1)) 

m r 

= - + ••• + ^- + ^ (5.(1) -re„(0)). (20) 

TO 1 j,a 

When a = 0, 

ea(l) - '^^.(O) = - - 2 + i (g„(0) - r ec(m - 1)), 
TO r 

and since ^^(l) = ^q(to — 1), we obtain 

gM) = --^ + ^-^^5cM 

TO 2 

aa(l)-rga(0) = l-l+ '^^+J ' g^{0) (21) 

TO 2 

From this point, the reader who wishes to verify details is encouraged to employ any standard 
computer algebra system. Equations (20) and (21) define a recurrence with an initial condition on 
differences of ^q, which is resolved by substituting ^/q(1) — r-C/c((0) from (21) into (20) and simplifying 
the geometric series. For a > this yields 

1 - r'' 1 / 1 
r'^~^m 1—r r" \m 2 



g.ia + l)-r g^{a) = — r-^— + - ( ZT - 1 + ^ ^"W ) • (22) 



Denoting the right-hand side of (22) by Da{a), for a > we have 
ea(a) = rg^ia-l) + D,,{a-l) 

= r''g,(0) + r"-iD„(0) + r"-2^„(l) + --- + r°D„(a-l). (23) 
A careful but straightforward summing of geometric series in (23) yields, for a > 0, 

1 l + r^"" 

g.{a) = -g^{o)- ^ 
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m r°"2 1— r\l — 1— ?' 1 + ?' 2r / 

Now using (24), we set Qa{^) = Gai'm- — 1) and solve for Qa{Q), obtaining 

C m - ^ ^ ^ 2r(l + r-) 

W - ^ _ 1)2 + (^2 _ _ 1) ' 

which together with (24) and simphfication yields 

2 _ ^-a/2-J ^^m/2-a/2 _ ^-(m/2-a/2) ) 



_2 2 (7-™/2 _|_ J,— JTi/2) 



T, (r + — 2) (r — r^^) (r™/^ — 7.-m/2) 

2 ^^a/2 _ ^-a/2) ^y,m/2-a/2 _ ^- (m/2-a/2) ~J 

_2 2(r™/^~l'^^^l + r^^/^+l^^^^l) 



m(r + r-i-2) _ j,-i-)(^m/2 _ ^-m/2) 

Substituting jy — x| for a gives the desired formula for Qa. □ 

By the definition of a and r, we may use the substitution r = e'^ to rewrite (r^ + r~'^) /2 = cos z0 
and (r^ — r~^)/2i = sinz0. Together with the definition of the Chebyshev polynomials of the first 
and second kinds, r„ and Un, respectively; i.e., 

Tn{x) cos nO and 
where x = cos0, wc obtain the following corollary to Theorem 8. 

Corollary 5. Let m > 3. For complex a ^ and < x,y < m — 1, the generalized Green's function 
Qa for the cycle Cm with vertices 0, 1, . . . , to — 1 satisfies 



ma a(2 + a)/7„/2-i (1 + Q^) 
where T and U are the Chebyshev polynomials of the first and second kinds, respectively. 

Proof of Theorem 7: The theorem follows by applying the integral formula for products of graphs 
without boundary in Corollary 3 to the torus, where T = Cm, T' = C„, the (/)"s are the orthonormal 
basis described in Lemma 1, Q and Q' are given by Theorem 1, and Qa is given by Corollary 5. □ 

Combining Theorem 7 with (3) using the orthonormal eigensystem of Lemma 1 for both Cm and 
Cn yields the following nontrivial identity. 
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Corollary 6. Let m, n > 3; < x, y < m — 1; and < x' , y' < n — \. Then 

1 exp ((27rij7TO)(y — x)) exp ((27rifc/n)(j/' — a;')) 

^ (1 - cos (2^j7m)/2 - cos (27rfc/n)/2) 



cxp {{2nik/n)[y' — a;')) 



fc=i 



1 



m(l — cos (27rfc/n)) 
r„/2-|y-x|(2 - cos (27rfc/7i)) 



(1 — cos (27rfc/n))(3 — cos {2iTk / n))U ^ /2-i{'^ ^ cos {2'Kk/n)) 

(m + l)(m — 1) Ij/ — xj (y — x)^ 

■ 5 1 

6mn n mn 



(n+l)(7i-l) \y'-x'\ [y'-x'f 



3mn m mn 

where T and U are the Chebyshev polynomials of the first and second kinds, respectively. 



The Laplacian of C,„ has a 1-dimcnsional eigenspace corresponding to eigenvalue 0, and a second 
1-dimensional eigenspace corresponding to eigenvalue 2 iff Cm is bipartite (when m is even). Oth- 
erwise, all eigenspaces are 2-dimensional, since Xj = Xm-j for all 1 < j < to — 1. This means that 
Corollary 6 is only one of a class of identities constructed by choosing orthonormal eigensystems for 
Cm and Cn for the left-hand side, and a possibly distinct orthonormal eigensystem for C„ on the 
right-hand side. 



5.2 The t-dimensional torus, Cmi x C^a x ■ • ■ x Cmt 

The bottleneck in computing the normalized Green's function for the t-torus via Corollary 3 is the 
lack of a formula for the generalized Green's function for any t > 2. Thus the decomposition into a 
cycle, for which Qa is given by Corollary 5, and a — l)-torus is required. 

Before giving the normalized Green's function of Cmi x • • • x Cmt > present the information 
on the components still needed. Choose and label the eigensystem of each Cm^ by 

{(A(:\0;.:^):O<j, <m,-l}, 

where m^ > 3, = Ag'''' < A^'''' < •■• < a|,j^^_]^, and the vectors {0^*'' : < js < — 1} are 
orthonormal. The eigenvalues of Cm^ x • ■ • x Cmt ^-^e averages of the eigenvalues of the factors 
Cm^, and the corresponding eigenvectors are products of the eigenvectors of the factors. This is 
summarized in the next well-known lemma, whose proof is a straightforward induction on (17). 

Lemma 2. The eigenvalues of Cm2 x ■ • ■ x Cmt '^'"^ 

where < js < m^ — 1 for all 2 < s < t, with corresponding eigenvectors 

ds)i 

s=2 
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For the following theorem, let Q be the normalized Green's function for Cmi from Theorem 1, 
and Q' the normalized Green's function for Cma x • ■ • C,„f 

Theorem 9. Let t >2. Let < Xj^, yj^ < rus — 1 where nis > 3 for 1 < s <t. The t- dimensional 
torus C'mi X • • • X C,„j has normalized Green's function 

G{{xi,...,xt),{yi,...,yt)) = t ^ ^k{x2, ■ ■ ■ ,xt)^K{y2, ■ ■ ■ ,yt)gAK{xi,yi) 

/^#(0,...,0) 

G'((x2, . . . (y2, ■ • ■ ,2/t)) + G{xi,yi), (26) 



[t — 1 jmi m2 • ■ ■ mt 

where K ranges over all indices (j2, ■ ■ ■ ,jt) 7^ (0, ■ • ■ ,0), $ and A are defined in Lemma 2, Qa is 
given by Corollary 5, and G' is the normalized Green's function for 

X ■ ■ ■ X Cm* ■ 

Proof. The proof proceeds by using F = Cmi and F' = Cm2 x • • • x Cmt in Corollary 4. The degree 
of F is d = 2, and the degree of F' is d' = 2(t - 1). The result follows. □ 

Although G' in Theorem 9 may not already be known, it can be computed inductively from 
repeated applications of the theorem. Determination of Qa for any small product Ci x • • • x Cf 
would allow the reduction in the number of applications required by a factor of f . The following 
formula for the 3-torus is obtained by two applications of Theorem 9, first taking the product of 
r = Cm with F' = Cm X Cm, and then the product of F = C,„ with F' ~ C„ 



^m- 



Corollary 7. For < j/i, X2, y2, 2^3, 2/3 < m — 1 where m > 3, the 3-dimensional torus Cm. x 
Cm X Cm has normalized Green's function 

^''((.Ti,.T2,a;3), (2/1,2/2,2/3)) = X! [exp ((27rzj7m)(2/2 - X2)) 

exp ((2nik/m){y3 - a:;3))^(l-cos(27ri/m)/2-cos(27rfc/m)/2)(a;i,2/l) 

3 ™ 

+ —^CXp{{2TTik/m){y3 - X3))5(i_eos (27rj/m)) (2^2, 2/2) 

3 

+ z;2 ( '"^ ' X^'^ ^' - l2/3 - 2:3! + 



m 
3 

3 



(m H 


- 1)(to- 


1) 




6m 




(m H 


~ 1)(to- 


1) 




6m 




(m H 


~ 1)(to- 


1) 


6m 



12/2 - 2:21 

|yi - 1 



(ys 






m 


(y2 






m 


(yi 


-x^f 


m 



where Qa is given in Corollary 5. 



These compact formulas for Green's functions of tori offer fast alternatives to computing pseudo- 
inverses of their Laplacians directly. This increase in speed, essentially due to the symmetry of the 
torus, is reflected in O(logn) computational complexity of the Chebyshev polynomials r„ and J7„. 
Various algorithms for computing r„ and Un are given in [10], and a more theoretical treatment of 
types of polynomials computable in O(logn) appears in [9]. The following corollary to Theorem 9 
is significant because the Laplacian (£, L, or A) of the torus on n vertices has rank n — 1, and so 
computing its pseudo-inverse provides along the way the inverse of an (n — 1) x (n — 1) matrix. 
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Corollary 8. Matrix pseudo-inversion of the Laplacian of the t-dimensional torus with n vertices 
via its Green's function is 0(t • n^^^/* logn), provided that the matrix itself is not completely recon- 
structed. 

Proof. We assume the t-dimensional torus is Cmi x • • ■ x C^t , where nig > 3 for 1 < s < t and 
nl=i ™s = fT-- It suffices to compute one row of the pseudo-inverse of the normalized Laplacian in 
order to know the entire inverse, due to the translational symmetry of the torus; i.e., since 

G((xi, . . .,xt), (yi, . . .,yt)) = G((0, ... ,0), {\yi - xi\, ...,\yt- xt\)). 

(In fact, because \ys — Xs\ can be replaced by to — \yg — Xs\ (mod m) without changing the value 
of G, only n!=ir"^s/2l of these entries must actually be computed.) Without loss of generality, 
"m-i > ■■■ > ^t- Then the summation term on the right-hand side of (26) has at most 
summands, which can each be computed in 0(t logn) time. The second and third terms can also be 
computed in 0(tn"'^~^/* logn) time, and we must compute n terms total to know all entries of the 
pseudo-inverse of £. □ 

For example, the time complexity of computing the Green's function for the 3-torus with n 
vertices is 0(n^/^ logn). Such a quick pseudo-inversion formula is surprising, since matrix inversion 
in general has the same complexity as matrix multiplication (see [2]). Matrix multiplication is known 
to be 0{n'^), where 2 < u; < 2.376 (see [5]). Surprisingly, for large n we can compute all of the 
values for the pseudo-inverse of the Laplacian for the torus faster than we can write down the whole 
matrix. Of course, requiring the presentation of the whole matrix rather than just the first row 
increases the complexity to 0(i logn -I- n^). 



6 Hitting times from Green's functions 

The equivalence of Green's functions to the fundamental matrix Z of (5) under similarity transfor- 
mation allows many quantities in random walks to be computed using Q. The hitting time Q{x, y) 
of a simple random walk starting at vertex x with target vertex y is the expected number of steps 
to reach vertex y for the first time by starting at x and at each step moving to any neighbor of x 
with equal probability. In [3], Chung and Yao show the following relationship between Q and Q. 

Theorem 10 (Chung, Yao). The hitting time Q{x,y) satisfies 

Qix,y) = ^-^giy,y)~^gix,y). 

We immediately have a computational formula for the hitting time whenever Q is known. For 
instance, whenever F is regular with n vertices, 

Q{x,y)^n{Q{y,y)-g{x,y)). (27) 

Figure 1 plots hitting time in the case of C4g x C4g, using Theorem 7 and (27). The domain consists 
of the vertices of the torus laid out in a square grid with periodic boundary, making opposite ends 
adjacent. The vertical axis plots the hitting time of a random walk initiated at (0, 0) with target 
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(x, y), achieving a minimum of at (0, 0) and leveling off at just above 6000 steps to reach vertices 
farthest from the start of the walk. If (27) is computed using techniques such as Corollary 5, 
the hitting time expression will involve orthogonal polynomials. Aldous and Fill claim in [1] that 
orthogonal polynomials appear whenever the graph has sufficient symmetry, but the dependence 
remains largely unstudied. 




Figure 1: Graph of the hitting time of a random walk on the torus C49 x C49, laid out as a square 
with periodic boundary 
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